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Abstract 

The effective theory of adiabatic fluctuations around arbitrary Friedmann-Robertson- Walker 
backgrounds — both expanding and contracting — allows for more than one way to obtain scale- 
invariant two-point correlations. However, as we show in this paper, it is challenging to produce 
scale-invariant fluctuations that are weakly coupled over the range of wavelengths accessible to 
cosmological observations. In particular, requiring the background to be a dynamical attractor, 
the curvature fluctuations are scale-invariant and weakly coupled for at least 10 e-folds only 
if the background is close to de Sitter space. In this case, the time-translation invariance of 
the background guarantees time-independent n-point functions. For non-attractor solutions, any 
predictions depend on assumptions about the evolution of the background even when the per- 
turbations are outside of the horizon. For the simplest such scenario we identify the regions of 
the parameter space that avoid both classical and quantum mechanical strong coupling problems. 
Finally, we present extensions of our results to backgrounds in which higher-derivative terms play 
a significant role. 
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1 Introduction 



A key feature of the primordial seed fluctuations observed in the cosmic microwave background 
(CMB) and the large-scale structure (LSS) is their scale-invariance [IJ [2]. Furthermore, the 
observed near-Gaussianity of the data requires the fluctuations to be weakly coupled^ over a large 
range of scales. In this paper we show that if the universe was dominated by a single degree of 
freedom at the time when the fluctuations were created, the above two facts together strongly 
constrain the background spacetime at that time. 

Our basic point is very simple: the scale-invariance of two-point correlations does not guaran- 
tee scale-invariance of interactions. Observations require scale-invariance of the power spectrum 
of curvature perturbations over at least AiV ~ 10 e-folds, from CMB scales (~ 10 4 Mpc) to 
galactic scales (~ IMpc). These fluctuations are created while the conformal time r evolves by 
a factor of e~ AN . In quasi-de Sitter backgrounds, this exponential change in time appears only 
in the time-evolution of the scale factor a(r) = — (Ht) , which is not a physical observable. 
All observable couplings during inflation are nearly time- independent. This is a consequence of 
the time-translation invariance of the background. The scale-invariance of the fluctuations pro- 
duced by inflation therefore applies to all n-point functions. In contrast, in non-de Sitter back- 
grounds with a single fluctuating degree of freedom^] (e.g. [6J [7]), it is hard to achieve both scale- 
invariant two-point correlations and time-independent interactions simultaneously. This is the 
case because the n-point functions of adiabatic (or 'single-clock') fluctuations around Friedmann- 
Robertson- Walker (FRW) backgrounds aren't independent, but are related by the symmetries 
of the background [U [9]: time translations are spontaneously broken in FRW backgrounds and 
the symmetry is non-linearly realized by a Goldstone boson. This non-linear realization of time 
translations forces unavoidable relations between the quadratic terms and the higher-order terms 
in the Lagrangian. The interactions are particularly constrained when the Goldstone boson of 
time translations is the only relevant degree of freedom. It is precisely these constraints that will 
allow us to make statements about the interaction Lagrangian in any cosmology that produces 
a scale-invariant two-point function. We find that any time-dependence that is not in the scale 
factor typically leads to an exponential growth of interactions, so that one has to be worried 
that the fluctuations become strongly coupled at some point during the evolution. At this point 
perturbative control of the theory would be lost. 

In Section [2] we review the conditions for scale-invariant two-point correlations in the effec- 
tive theory of adiabatic fluctuations. At lowest order in derivatives, we identify two classes of 
backgrounds that allow for scale-invariant two-point functions for the primordial curvature per- 
turbation C- These two possibilities are distinguished by whether £ is constant on superhorizon 
scales (Case I) or grows in time by an exponentially large amount (Case II). These two possibili- 
ties correspond to the background being a dynamical attractor or not. Within these two classes 
of solutions there are important subclasses depending on whether the background is expanding 
or contracting, has a strongly time-dependent equation of state or a strongly time-dependent 

ll Weak coupling' implies that the action for the fluctuations has a well-defined expansion in which the leading 
terms of order n + 1 in the fluctuations are smaller than the terms of order n. We shall be more precise about this 
in 931 




11 of our constraints can be evaded in multi-field non-de Sitter cosmologies such as [3] 111 E] . 
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speed of sound. In Section [3] we explain that higher-order interactions of £ and the associated 
higher-point correlations severely limit the possibilities. We present the leading cubic interac- 
tions, identify the time-dependent 'coupling constants', and discuss when they lead to a strong 
coupling problem. In Section [4] we study the attractor solutions in detail. For the case of a 
constant speed of sound we find exact solutions which allow us to study FRW backgrounds in 
full generality. We show that only near-de Sitter backgrounds, with intrinsically small time- 
variations of all quantities, avoid the strong coupling problem. We explain why our conclusions 
are not affected significantly when the speed of sound is allowed to vary in time. We present our 
conclusions in Section [H 

Four appendices contain further technical details: 

In Appendix |A| we derive the quadratic action for the curvature perturbation £ in the frame- 
work of the effective theory of adiabatic fluctuations [8j |9|. We show that, at leading order in a 
derivative expansion, the dynamics is characterized by two free functions: the Hubble rate H(t) 
and the sound speed c s (t). In the limit in which M^H becomes smaller than other parameters 
in the Lagrangian, higher-derivative terms can be the leading effects. We discuss this limit in 
Appendix [B} In Appendix [C] we turn our attention to the non-attractor cases. These theories 
are significantly less predictive and require a variety of additional assumptions about the evo- 
lution both before and after the scale-invariant modes are produced. We will show that even 
granting the most optimistic such assumptions the parameter space is significantly constrained 
by the weak-coupling requirement. Finally, in Appendix [D] we compute the full bispectra for the 
non-attractor cases of theories with non-trivial speed of sound, following earlier work by Khoury 
and Piazza |10j . 



2 Scale-Invariance of the Two-Point Function 

Although originally introduced to study the possibility of violating the null energy condition [8] 
or to describe adiabatic fluctuations around inflationary backgrounds |9], the effective theory of 
adiabatic fluctuations applies to arbitrary FRW backgrounds. It is therefore a promising starting 
point to discuss the properties of fluctuations around a wide class of background spacetimes, both 
expanding and contracting. Using this framework we show in Appendix [A] that the quadratic 
action for curvature perturbations £ (at leading order in a derivative expansion of metric pertur- 
bations in comoving gauge) depends only on two free functions of time: the scale factor a(t) and 
the speed of sound c s (t). In addition, we have the following derived quantities 



H 

IP 



V 



e 



He, 



(2.1) 



where dots are derivatives with respect to time t and H = dt In a is the Hubble parameter. During 



inflation the parameters in (2.1) are all much less than unity and the speed of sound is nearly 



constant, but here were want to consider a broader spectrum of possibilities. 

Our starting point is the quadratic action for curvature fluctuations (see Appendix |A|) 



S 2 



3 

Mpj / dtd 3 x^ 



c 2 



r(^C) 2 



(2.2) 
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Rescaling time, 



this can be written as 



where 



dy = c s dr 



-dt , 



S 2 = M 2 l Jd*x q 2 [(C') 2 



a 2 e 



(2.3) 
(2.4) 

(2.5) 



The prime in (2.4) indicates a derivative with respect to y and d 4 x = dyd 3 x. The time variable 
y is convenient because a mode of comoving wavenumber k crosses the (sound) horizon when 
|A;y*| = 1. We are interested in the regime where modes cross the shrinking horizon \y\, so 
that we concentrate on negative y. The behavior of Ck(y) changes from oscillating modes to 
monotonically growing and decaying modes at y*. The prefactor q 2 in (2.4) has to be positive: a 
negative q 2 implies a wrong-sign kinetic term for £ and hence a ghost instability. Notice that a 
negative q 2 implies a violation of the null energy condition (H > 0). This can be achieved only 
using higher-derivative terms in the effective action [8]. 

The time-dependence of the prefactor q(y) determines whether the fluctuations have scale- 
invariant two-point correlations. To see this, let us define the canonically-normalized variable 
v = \^2q(. After a few integrations by parts the action becomes 



S-2 



M 2 pl 



d*x 



./\2 



(d t v) 2 + 



a" 

" „.2 



This corresponds to the following equation of motion in Fourier space 



where 



vl+ k 
n(n — 1) 

y 2 



q 



Vk = 0, 



(2.6) 



(2.7) 



if 



oc y 



(2- 



For the power law ansatz (2.8), the differential equation (2.7) has an exact solution in terms of 
Hankel functions of the first and second kind: 

v k {y) = [cxH^i-ky) + c 2 Hl 2 \-ky)] , 



where v 2 — — = n{n — 1) . (2.9) 



Demanding that the solution reduces to the standard Minkowski limit on small scales (or at early 
times) , 



lim v k (y) 



-iky 



fixes C2 = and c\ 



|. The power spectrum on superhorizon scales then is 



lim fc 3 |ufc| 
ky— >0 



1 



(-ky) 



3-2u 



(2.10) 



(2.11) 
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where we used lim^^o Hjp ~ (—ky) u . This describes scale-invariant fluctuations if and only if 



v-- = 2 = n(n-\) 
4 ; 



n = {-l, 2} . 



(2.12) 



We hence found two distinct cases for which the background q(y) leads to a scale-invariant 
spectrum of fluctuations [T0]:[^] 



Case I 



Case II 



q = Qi 



(2.13) 
(2.14) 



On superhorizon scales the growing mode of £ is constant in Case I and time-dependent in Case II: 



lim £ 

ky— >o 



lim — 

ky^O q 



1 

'3L)~ 



Case I 
Case II 



(2.15) 



This difference in the time-dependence of C, reflects the (in)stability of the background. In Case I, 
C approaches a constant on superhorizon scales and the background is a stable attractor solution; 
in the long- wavelength limit, k — )■ 0, classical solutions for £ are interpreted as homogeneous 
perturbations to the background scale factor: ds 2 = a 2 (r)(— dr 2 + e 2 ^ T ' x )dx 2 ). In contrast, the 
exponentially large growth of £ in Case II indicates that the background is unstable and not an 
attractor. This instability implies that Case II is not predictive unless additional assumptions 
are made about the evolution before and after the scale-invariant modes exit the horizon (see 
Appendix [Cl). 



Several examples of backgrounds satisfying (2.13) and (2.14) are known: for constant speed 
of sound c s = const, the condition for scale-invar iance, 



a 2 e 



a 



a"a 



!\2 



OC T 



(2.16) 



implies that either the scale factor a or the equation of state e (or a combination of both) 
vary rapidly with time. The case a oc — r _1 and e ~ const, of course corresponds to slow-roll 
inflation [T3], while e oc r~ 2 and a ~ const, corresponds to adiabatic ekpyrotic contraction (or 
expansion) with time- varying equation of state [6j [7] . In Section [4] we will find an exact solution 



to (2.16) which to our knowledge is new. This solution contains inflation and adiabatic ekpyrosis 



as special cases. Allowing for a non-trivial speed of sound, Khoury and Piazza [10J discussed the 
production of scale-invariant fluctuations in non-attractor backgrounds. Because of the many 
extra assumptions and model-building requirements, we find the non-attractor solutions much 
less compelling than the attractor solutions. In the main text we will therefore focus on the 
attractors cases, leaving a detailed treatment of the non-attractor cases to Appendix [Cj 



3 These solutions may also be found without assuming q" /q to be a power law, by requiring invariance of the 
action (2.4 1 under the rescaling y — > Ay, x — > Xx and £ — > A n £, with n some number. If the initial state does not 



break this symmetry — as is the case for the Bunch-Davies vacuum — then the two-point functions of fields with 
scaling dimension are scale-invariant. 
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3 The Strong Coupling Problem 



The main point of our paper is the following: although there are many different FRW back- 
grounds - both expanding and contracting - that lead to scale-invariant two-point correlations 
(see Section [4] and Appendix [C]) , in most cases instabilities are revealed when looking at higher- 
order correlations. In other words, scale-invariance of the two-point function does not guarantee 
scale-invariance of the interactions. This is the case because the non-linear realization of time- 
diffeormorphisms forces non-trivial connections between the quadratic and the cubic Lagrangian. 

The cubic Lagrangian contains among others the following terms [15] (see Appendix [Pj) 



C 3 c 



a 3 e 



(dC) 2 1 

(e - 3(1 - c 2 s )) CC 2 + (e - 2e s + 1 - c*) C~^- + ^{r 1 + Hr ] - 2He s ) ( 2 ( 



(3.1) 

where S3 = J dtd^x £3. As a simple measure of our strong coupling concern we use the ratio of 
cubic to quadratic Lagrangian: 



X 



£2 



0({l,e, V ,e s }) 



c 



< 1 



(3.2) 



When this ratio (or more generally the ratio becomes larger than one, the theory is 

strongly coupled. There are two different ways in which having C- in t/ C2 larger than one is dan- 
gerous. The first is when this condition happens at horizon crossing. In this case the theory is 
strongly coupled at the quantum level, in the sense that quantum loop corrections aren't sup- 
pressed [9J [16] . This is the regime which is usually referred to as strong coupling. In the absence 
of a UV-completion or a proof of strongly coupled unitarization, unitarity is lost and the model 
is not under theoretical control. For attractor solutions (Case I), this is the only condition that 
we will require. 

For non-attractor solutions (Case II), the condition (3.2) can also be violated after horizon 
crossing. In this case, C- m t > £2 does not represent a violation of unitarity or an actual quantum 
mechanical strong coupling. It simply signals the presence of classical non-linearities. In principle, 
these can be treated by solving the non-linear equations for the modes. Hence, strictly speaking, 
the model is still under theoretical control. However, the predictions that are made by studying 
only the quadratic Lagrangian, including the prediction of scale invariance, can't be trusted 
anymore. Moreover, if such a regime happens for the observable modes, then the fluctuations can 
be strongly non-Gaussian and in conflict with observations]^] We will need to use this classical 
strong coupling condition only when studying the non-attractor solutions (Case II) in Appendix [C] 

In Case I, £ freezes at horizon-crossing, so that its magnitude at that time is fixed by ob- 
servations, ~ 10 -5 . Similarly, the speed of sound at horizon-crossing can't be too small, 



> c 



1/2 



10- 



to avoid strong coupling from interactions proportional to c s ( 



straint X < 1 then becomes a constraint on the size of the 'couplings' e, r] and e s 



). The con- 
2. If these 



4 Recall that to be consistent with present microwave background and large-scale structure observations the 
background has to allow for at least of order AN w 10 e-folds of scale-invariant and weakly-coupled fluctuations. 

5 Observational constraints on non-Gaussianity [12] of course imply an even stronger constraint on the minimal 
speed of sound, c 3 ,+ > 0.01. 
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couplings are time-dependent one has to worry that the theory becomes strongly coupled at some 
time in the regime of cosmological interest. We discuss this in detail in Section |4j 

In Case II, the analysis is significantly complicated by the fact that ( evolves outside of the 
horizon. Its value at horizon-crossing, can be exponentially smaller than its value at the 
end of the scaling regime, Cend- Whether the fluctuations are weakly or strongly coupled at 
horizon-crossing then depends on the time-evolution of the speed of sound. Furthermore, the 
super horizon evolution of £ can lead to additional classical non-linearities. In Appendix [C] we 
discuss the combined constraints from non-linearities at horizon crossing and the subsequent 
classical evolution. 



4 Scale-Invariance for Dynamical Attractors 

Only in Case I is the background a dynamical attractor. This means that curvature fluctuations 
freeze at horizon crossing and the predictions for observations are independent of the subsequent 
evolution. We will first consider the limit of constant sound speed, c s ~ const., and then explain 
why a time-dependent sound speed does not qualitatively change our conclusions. 



4.1 Emden- Fowler Equation 



As we have shown above, theories with constant sound speed produce scale-invariant fluctuations 
if the scale factor a(r) satisfies the following non-linear second-order differential equation, 



a 2 e 



a 



(a') 2 



,2 <i 



" r 2 



(4.1) 



where ' denotes derivatives with respect to conformal time r. We assume that (4.1) is valid in a 
finite time-interval between n and r en d- During this 'scaling regime' the background is consistent 
with scale-invariant two-point correlations. It is convenient to normalize the scale factor such 
that it is unity at the beginning of the scaling regime, 



1 



(4.2) 



This can always be achieved by a conformal rescaling of the spatial three-metric. It will also be 
helpful to define a dimensionless time variable 



HiT, 



where Hi = H(ri). By a change of variables, b = 1/a, equation (4.1) becomes 

a 



b" 



:(b')% 



(4.3) 



(4.4) 



where ' now denotes differentiation with respect to z and we have defined the constant a = c s tizf. 
Equation (4.4) is a special case of the generalized Emden-Fowler equation |ITJ, whose exact 
solution we present below. The 'slow-roll' parameters derived from b(z) are 

b"b _ a 2 (hie)' _ 2 b 



2 . 



(4.5) 
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By definition, cf. (4.2) and (4.3), the function b{z) satisfes the following boundary conditions: 



b(zi) = 1 and b (zi) = —1 



(4.6) 



The solutions to (4.4) are classified by the two constants Z{ and ej (or a).^ While = e(zi) 
characterizes the equation of state at the start of the scaling regime, the variable zi = H{Ti is the 
initial ratio between the freeze-out scale, |tj|, and the (comoving) Hubble scale, (aj-ffj) -1 , 



Zi 



(aiHi)' 



(4.7) 



In quasi-de Sitter spacetimes the two scales coincide, Zi = — 1 — e% ~ —1, while more generally 
(i.e. for time-dependent e) the scales can be quite different. We consider the branch of solutions 
with T{ < Recalling that |r| decreases by at least a factor of e~ 10 during the scaling regime 
implies \z en( ±\ < e _10 |zj|. This exponential change of conformal time during the scaling phase will 
feature prominently in the strong coupling problem. 



4.2 Approximate Solutions 



To develop some intuition for the space of solutions to (4.4), we will first derive approximate 



solutions as an expansion in small a. This reproduces three known limiting cases: slow-roll 
inflation [13] . slow contraction with a rapidly changing equation of state [6], and slow expansion 



with rapidly changing equation of state [7]. Then we will show that (4.4) in fact permits an exact 



solution (£ |4.3[ ). This will allow us to study FRW backgrounds in full generality (^ ]4.4[ ). 

An iterative solution to (4.4) as an expansion in powers of a (rh may be written as follows: 



where the functions b 



iff 

b { 2) 

&/ (3) 



(n) 

, 

Q 

a 
a 

^2 



K z ) = ^2 b (n)(z) 

0(a n ) satisfy 



n=0 



\0)) 0(0) 



( fe (o)) 2fe (i) + 2b {o) b {i) b (o) 
(6( )) 2 6(2) + 2 6(0)6(2)^(0) + (6(i)) 2 6( ) + 2 6 , ( 0) 6 / ( 1) 6 (1 ) 



(4.8) 

(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 



The expansion converges if 6(„ +1 ) <C 6( n ) and 6( n+1 ) -C biy The initial conditions (4.6) imply 
6(o)(^) = -6( )( z i) = 1 and b (n>i){zi) = 6'( n > 1) (z i ) = . (4.14) 



6 From now on we will restrict ourselves to the limit c s f» 1, as at this point it is clear that the solution for 
c s w const, can be trivi ally r ecovered from the solution for c s = 1 by substituting r — > c s r and a — > c 3 a. 
7 Note that equation (4.1 1 has a Z2 symmetry, z — > —z, relating the solutions for r, < and r, > 0. 



Strictly speaking our solutions will be expansions in a/z and will break down for small \z\ < a <C 1. At that 



point we will switch to the exact result of [ 4.3 
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To zeroth-order in a, the solution to (4.4) therefore is 



5 (0) 



l + Zi 



(4.15) 



We find the first-order correction bm by substituting (4.15) into (4.10), integrating (4.10) and 
imposing (4.14): 



a: 



(1 + Zi + z) In - + (1 + 2 Zl ) (l - -) 



Hence, the solution to first-order in a 



1 + (zi - z) - Eizj 



uzf is 



;i + Zi + z)ln- + (l + 2zi)(l 



(4.16) 



(4.17) 



Continuing this iterative procedure we could determine the solution to arbitrary order in a. 
However, all the interesting physics is already included at this order of the approximation. In 
particular, we will now show that (4.17) has two limits of particular interest: slow-roll inflation 
and slow contraction or expansion with rapidly changing equation of state. 

Slow-roll inflation. In quasi-de Sitter space the ratio of the freeze-out scale to the Hubble 
scale is 



In this case, the result (4.17) reduces to 

b = - 

Using 

H 

M = 



_l_ e . + 0(e 2) . 

-z[l - €i + a ln\z\ 
■■ -b' = 1 + e* In - 



this implies the following solution for the scale factor 



1 

Ifr 



(4.18) 

(4.19) 
(4.20) 

(4.21) 



This is the conformal time for a quasi-de Sitter space, a familiar result from the standard treat- 
ment of inflation |13|. 



We should be careful about the range of validity of the solution (4.19): our convergence 
criterium, <C fy ) an d &m ^ ^(o)> w ^ n break down for small z. In that limit the stronger 
condition is bn\ <C br \. Specifically, for \z\ <C \z%\, (4.15) and (4.16) become 



'(0) 



-e„: - z 



and 



&(i) 



(4.22) 



and we hence trust our approximate solution in the regime e« < \z\ < \zi\. In this regime the 
solution (4.19) should be a good approximation to the exact answer. 



Substituting (4.19) into (4.5), we infer the slow-roll parameters 



e = €i(l + 2 £i ]n\z\) + 0(e 



1 1 i 



(4.23) 
(4.24) 
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This result is consistent with the fact that at first-order in a slow-roll expansion the deviation 
from scale-invariance is given by n s — 1 = — 2e— r/. At first-order in ej, we indeed find a cancellation 
between e and 77: n s — 1 = C(e|). At second order in ej there is a small time-dependence of e and 
r\. By construction, this time-dependence is exactly such that scale-invariance is preserved. 



Finally, we recall the weak-coupling criterium (3.2) 



X = {e, V }C, < 1 



(4.25) 



From the time-independence of the couplings (4.23) and (4.24) we see that the slow-roll solution 

(4.26) 



is indeed weakly coupled for sufficiently long times 

X(z) <C 1 for Zi < Z < Zend 



Slow contraction. For positive z (i.e. negative H{) the result (4.17) describes contracting 



solutions. Furthermore, in the limit < Zi <C 1, the solution (4.17) becomes 



1 + Zi 



tiZ; 



In h 

Zi 



(4.27) 



As before, this is a good approximation to the exact solution as long as <C 6(0) and b'^ <C b'^ 
Now the stronger condition is b 1 ,^ < b', Q ^ . Since 



'(o)- 



bi 



(0) 



and 



(i) 



2 i 



we therefore trust the approximate solution (4.27) in the regime eizf < 



\Z\ < \Zi 



'slow- roll' parameter derived from (4.27) grows exponentially with time 



(4.28) 
The first 

(4.29) 



Essentially, since the scale factor is nearly constant in the solution (4.27), all the time-dependence 
of q 2 = a 2 e = e/b 2 oc 1/z 2 is provided by e rather than by a. This is the contracting adiabatic 
ekpyrotic phase first discussed by Khoury and Steinhardt [6]. 

Since the time coordinate z changes exponentially - by at least a factor of e~ 10 during the 
scaling regime - one has to worry about the exponential growth of e as a source of large interac- 
tions. Strong coupling can only be avoided if ej is exponentially small [6], in which case there is 
not enough time for it to blow up before the end of the scale-invariant phase. However, in this 
case the second slow-roll parameter typically still leads to trouble: 

(hie)' _ _1 
'(In 6)' ~ z ^ 



\m 



> 



1 . 



(4.30) 



Independent of the size of the coupling n is large and exponentially growing. In particular, 
after less than 10 e-folds the theory becomes strongly coupled, X ~ T)-kC* ~ 1- 

To show that this is a generic problem and not in any significant way tied to our expansion 



scheme, we consider the time-derivative of the e parameter in (4.5): 

V 



d_ 

dz 



(InVi) 



1 

-z + -b 



(4.31) 
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-H/Hi, this can be written as 

d 



Using b' 

r . 1 H 

During slow-roll inflation with H = —b/r > 0, a cancellation between the two terms on the 



(4.32) 



r.h.s. of (4.32) allows for a nearly time-independent e parameter. However, in a contracting 



spacetime with H < such a cancellation is impossible. Instead, we obtain the inequality [7] 

(4.33) 



f (Wi)>n 

ar t 



Hence, during a contracting phase e is forced to change exponentially, typically leading to a large 
77-coupling in the cubic action. 



Slow expansion. Since equation (4.4) has a symmetry under z —> —z, we expect a similar 
solution, still with a > 0, corresponding to slow expansion with rapidly changing e. Indeed, for 



negative z (i.e. positive Hi) the result (4.17) describes expanding solutions. Furthermore, in the 



limit \z{\ <C 1, the solution reduces to (4.27). This is the expanding adiabatic ekpyrotic phase 



discovered by [7J. Our above statements concerning the problem of large w still apply in this case. 
4.3 Exact Solutions 



We now show that the non-linear differential equation (4.4) in fact has an exact solution. Consider 
first the generalized Emden- Fowler equation 



b" = az n b m {b') 1 



(4.34) 



For {n, m, 1} = {—2, 1, 2} this is precisely of the form of (4.4). Remarkably, this special case has 
an exact solution 111 In: 



z(t) = f3 lCl [G{t)+c 2 ] , 
b(t) = (3 2Cl S.2t[G(t) + c 2 } + 



dG 
~dt 



where a = \fi\l3 2 2 and 



G(t) 



~ t2 dt 



Erf(t) 



(4.35) 
(4.36) 

(4.37) 



Digression: The following manipulations massage this result into a more user-friendly formp"| 



First, we may combine (4.35) and (4.36) into 

dz 



b = /3 



2tz + 



dt 



where 8 = = ± . — 



Furthermore, the first derivative of b{z) is 

db db 



b' 



dz dt 



dz 


-1 


dz 




= 2Bz 






A. 



(4.38) 



(4.39) 



Here, we correct a crucial typographical error in the result of [TT] 



10 Readers without an affinity for these basic algebraic manipulations may jump directly to [ 4.4 for a graphical 
representation of the solution. 
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where we used 4jf = — 2t4|. Given Zi and ej we can therefore map the initial conditions b(zi 



into the parameters c± and C2: We note from (4.391 that b'(zi) = —1 implies 

= -20* . 



dz 



The constraint b(zi) = 1 then gives 



U = j3 



1 



2/3 2 z, 



Equation (4.35) becomes 



= l- 7 [G(t)-G(ti)} , 



where 



7 



/?1C1 



2/3 e* 



To implement the solution numerically it is useful to define x = t — ti and hence obtain 



1 [G{t)-G{t l )] = 2[}e t 



t,2 dt' = 2/3 / er 2Ux '- x ' 2 dx' . 



Finally, we find 



Z(x) = Zi 



1-2/3 / t 




-ItiX -x 1 ' 



dx' 



-b'{zi) = 1 
(4.40) 

(4.41) 
(4.42) 
(4.43) 

(4.44) 

(4.45) 



b(x) = 2/3zi 



( U +x)^-(3e- 2t ^ x 2 



(4.46) 



These solutions have two input parameters: z% and j3 (or ej). Requiring \z\ to decrease with increasing |a;| 
imposes x > for (3 > and a; < for /3 < 0. In fact, recalling that ti is an odd function of f3, we discover 
that the solutions have the following symmetry 

Z{z t ,p}(%) = z {*i,-0}(~ x ) and b {z % ,p}{x) = b {Zi _p}{-x) . (4.47) 
Without loss of generality we can therefore focus on the /3 > branch of the solutions. 



From (4.45) and (4.46) it is straightforward to compute the slow-roll parameters using (4.5). With 

(4.48) 



b x = 2(3z we get 



V (x) = 2 z iJ_ e _ 2Ux _ x 2_ 2 

Z (-Z) 



4.4 Illustrations of the Strong Coupling Problem 



By construction, fluctuations living in the backgrounds (4.46) have scale-invariant two-point 



correlations. However, as we have argued above, this does not imply that higher-order correlations 



are truly scale- invariant. In this section we use the solution (4.46) to analyze the strong coupling 



problem of ^3] in a model-independent way. We note that our solutions are determined by two 
input parameters: z% — the ratio of the freeze-out scale to the comoving Hubble scale at the 
beginning of the scaling phase — and ej — the equation of state at that time. The parameter Z{ 
may also be viewed as a measure of the deviation of the solutions from the inflationary quasi-de 



Sitter backgrounds. For inflation we have %i - 
— 1 are very different from slow-roll inflation. 



-1 — ei + e|, while backgrounds with z± far from 
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In this section, we will scan over the possible values of the input parameters Zi and ej. As a 
diagnostic for the strong coupling problem we will then compute the rj parameter ./V e-folds after 
the beginning of the scaling regime: 



VN 



(4.49) 



where z = e z%. For logr/jy ^ 4 we lose perturbative control. Even before that, for logr/N > 2, 
we expect tension with constraints on non-Gaussianity. 

Slow-roll inflation. This is the most well-known case. We have perturbative control for an 
arbitrary amount of time. We have confirmed numerically that this regime is accurately described 



by the approximate solution (4.19) (within its expected range of validity) and the exact solution 



(4.45) and (4.46). For the exact solution to produce more than 60 e-folds may require Zi to be 



defined to second order in ej. 

Slow contraction. In contrast, in Figure[T]we show a canonical example of a slowly contracting 
spacetime with ej = 10 -8 and Z{ = 0.01. As expected, e and rj grow exponentially and within a 
small number of e-folds the theory becomes strongly coupled. At this point perturbative control 
is lost and, in the absence of a UV-completion or a proof of strongly coupled unitarization, the 
theory becomes unpredictive. 

log{e, T]} 




N 



Figure 1: Slow contraction: plot of log e and logry vs. N for Ci — 10 8 and Zi — 0.01. 



de Sitter vs. non-de Sitter. In Figure [2] we present a more general scan of the parameter 
space of solutions. For fixed ej = 10 -4 , 10~ 5 , 10 -6 , 10~ 7 we scan the parameter Zj. We plot the 
maximal value of Tj in the interval \z%^ -^end — ^ 

] . We find numerically that away from Zi — 1 
or | Zi | 1 it is hard to produce even just 10 e-folds of scale- invariant modes unless ej is very 
small (see also [7]). This is the reason that the points corresponding to e» = 10~ 4 and 10~ 5 don't 
extend to larger values of \zi\. Only for the range of z; t shown does the background produce at 
least 10 e-folds of modes. For z% ~ — 1 or \zi\ <C 1 the background allows more than 10 e-folds 
even without exponentially small ej. 
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log 



Vio 




-2 



-4 - 

Figure 2: The vaiue of the rj parameter 10 e- folds after the beginning of the scaling phase illustrates the 
strong coupling problem arising in backgrounds that deviate too much from quasi-de Sitter (zt — — 1 — e^J. 



The plot includes all three phases mentioned above: i) slow-roll inflation {z\ ~ — 1 — + e?), 
ii) slow-contraction (0 < zi -C 1) and in) slow-expansion (—1 <C Zi < 0). We see that all 
backgrounds except quasi-de Sitter suffer from a strong coupling problem, making it highly 
unlikely that more than 10 e-folds of scale-invariant modes can be realized in a controlled way. If 
one is satisfied with just producing the modes observed in the CMB, one may consider creating 
~ 5 e-folds of scale-invariant modes in a non-de Sitter evolution. To decide whether these modes 
would be consistent with observational constraints on non-Gaussianity would require a more 
careful computation of the three-point function including the effects of other interactions that 
could partially cancel the large contribution from the interaction proportional to r/. Why those 
few e-folds would happen to be the ones that we have observational access to would require an 
explanation. 

4.5 Generalizations 

So far our analysis was limited to actions with trivial propagation speed for the cosmological 
curvature perturbations, c s = const. In this subsection, we offer a few comments on extensions 
of our results to theories with non-trivial speed of sound and to cases in which high-derivative 
terms play an important role in the dynamics. 

Varying speed of sound. For attractor solutions, the requirement of weak coupling at horizon 
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crossing implies a lower limit on the sound speed 



C 2 



C 



< i 



1/2 



10" 



(4.50) 



The observational upper limit on the level of non-Gaussianities in the CMB strengthens this 
bound by one or two orders of magnitudes (12j f^~| Moreover, in order to avoid superluminality 
the sound speed is bounded from above, c s < 1. Hence, the sound speed can vary at most by one 
or two orders of magnitude over the entire range of the scaling regime from yi to y en d = e~ AN Ui- 
The allowed time- variation of the speed of sound is therefore too constrained to have a significant 
effect on g 2 oc l/c s . In particular, we do not believe that allowing the speed of sound to be 
time-dependent will give qualitatively new, non-de Sitter solutions. 

If we are looking for a scale-invariant solution that can last at least 60 e-foldings, then the 
interactions proportional to e and r\ will force us to have e approximately constant. In this regime, 
the requirement for scale invariance becomes: 



a 2 e 



c s ,i V y 



c s oc (ayf 



(4.51) 



This allows us to derive a(y) as a function of y: 



H 
IP 



a"a 



{a'f 



2 a 
ya' 



i 



Plugging this into (4.50), we get 



C 2 



4e 
1+e 



c 



(4.52) 



(4.53) 



To avoid that the interactions grow catastrophically as y — > 0, we require < e C 1 or e < 0. 
The second option makes the fluctuations ghost-like, so we are left with inflation as the only 
weakly-coupled solution. 

Higher-derivative terms. In Appendix [A] we derive the most general quadratic action for 
curvature fluctuations ( in the effective field theory approach of |9] . In unitary gauge there 
are no matter fluctuations, but only metric fluctuations, so we define the effective action by 
writing down all operators that are functions of the metric fluctuations and invariant under 
time-dependent spatial diffeomorphisms. The two most important elements appearing in this 
construction are the metric perturbation 5g 00 and the extrinsic curvature perturbation 5Kij. In 
Appendix [A| we use these geometrical quantities to define the most general action with unbroken 
spatial diffeomeophisms. As usual in effective field theories, this can be done in a low-energy 
expansion of the fields and their derivatives: 5g 00 is scalar with zero derivatives acting on it, 
while 5Kij is a one-derivative object. In many situations the terms involving 5g 00 therefore 
dominate the dynamics. This leads to a quadratic action for £ with non-trivial sound speed 
and is hence captured by the analysis in this section. However, in some limits, identified by a 



As we stressed before, having a small, constant speed of sound is not much of an issue once the constraint from 
non-Gaussianities is satisfied. 
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particular symmetry structure of the Lagrangian for the fluctuations, the terms involving the 
extrinsic curvature 5Kij can in fact become the dominant effects. For instance, in the limit in 
which Mpj-ff is smaller than a combination of the coefficients of the extrinsic curvature terms 
and the Hubble parameter H (see Appendix [B|) , the extrinsic curvature terms can lead to a de 
Sitter background with significant higher-derivative terms in the action for £. These are the ghost 
inflation [19] regime and the solutions identified as 'near-de-Sitter limit' in [9] and |12j . 

Can there be non-de Sitter backgrounds in which higher-derivative terms allow weakly-coupled 
scale-invariant fluctuations? In Appendix [B] we address this question. Although we don't give 
a comprehensive analysis of the range of possibilities implied by the rather complicated actions 
for £, we provide an extensive discussion of the most interesting limits of the parameter space. 
We reproduced the known near-de Sitter solutions, but found no non-de Sitter backgrounds with 
weakly-coupled scale-invariant fluctuations. 

5 Summary 

Cosmic microwave background (CMB) and large-scale structure (LSS) observations constrain 
the primordial seed fluctuations to be i) nearly scale-invariant and it) approximately Gaussian 
(or weakly coupled). In this paper we asked what these two observational facts together teach 
us about the cosmological background at the time when these fluctuations exited the horizon. 
We considered the most general effective theory of adiabatic fluctuations around completely 
arbitrary FRW backgrounds allowing both for expansion and contraction. We assumed that 
the observed cosmological perturbations arose from quantum mechanical fluctuations of a single 
degree of freedom. For theories with constant sound speed and requiring the background to be 
a dynamical attractor, we reduced the requirement of scale-invariant two-point correlations to a 
non-linear differential equation for the scale factor a(r), 



This equation has an exact solution — a special case of solutions to the generalized Emden-Fowler 
equation |llj . Using this solution we reproduced two known limits: slow-roll inflation [13J (a oc 
t -1 ) and ekpyrotic contraction (or expansion) with rapidly changing equation of state [61 [7] 



The non-linear realization of time diffeomorphisms in the effective theory of the fluctuations 
forces specific relationships between the coefficients of the quadratic Lagrangian and the in- 
teraction Lagrangian. This drastically limits the number of viable models that can produce a 
scale-invariant two-point function while staying weakly coupled for a sufficiently long time. In 
fact, we showed that only inflation leads to nearly time-independent couplings of all higher-order 
interactions such as e, r\ and c~ 2 . This implies that only inflationary spacetimes allow the fluctu- 
ations to stay weakly coupled over the range of scale relevant to cosmological observations. For 
all non-de Sitter backgrounds we identified a strong coupling problem that limits the predictivity 
of these solutions. Within 10 e-folds after the beginning of the scaling phase the perturbative 
expansion of the theory breaks down. This makes it challenging to produce the modes observed 
in the CMB and the LSS in a consistent non-de Sitter background. 
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Figure 3: The value of the r\ parameter 10 e-folds after the beginning of the scaling phase illustrates the 
strong coupling problem arising in backgrounds that deviate too much from quasi-de Sitter (z^ = — 1 — (4). 

Dropping the requirement that the background should be an attractor, we found non-de 
Sitter solutions with rapidly time- varying speed of sound (following previous work by Khoury and 
Piazza [TO]). In Appendix [C] we identified the regime of parameter space in which those theories 
evade both quantum mechanical and classical strong coupling problems. Since these backgrounds 
are not attractors, the curvature perturbations Q evolve after horizon exit and predictions depend 
on assumptions about the physics both before and after the regime during which the scale- 
invariant modes exit the horizon. We described the basic model-building requirements for non- 
attractor non-de Sitter backgrounds with weakly-coupled scale-invariant fluctuations. Whether 
these physical elements can be realized in a coherent theoretical framework remains an open 
question. In contrast, it is quite remarkable how the time-translation invariance of quasi-de 
Sitter backgrounds without any additional physical ingredients solves the horizon and flatness 
problems, while at the same time allowing for scale-invariant n-point functions. 

Note added. While this paper was being completed, Ref. |7j appeared which has some overlap 
with our Section 01 
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A Quadratic Action for £ 



In this appendix we derive the second-order Lagrangian for scalar curvature fluctuations £ around 
arbitrary FRW backgrounds, including the mixing of the scalar with gravity [8, 9j. We start in 
unitary gauge, in which there are no matter fluctuations, but only metric fluctuations. The 
most general effective action is then constructed by writing down all operators that are functions 
of the metric fluctuations and invariant under time-dependent spatial diffeomorphisms. Time 
diffeomorphisms are broken by the time-dependence of the background. For this purpose it is 
particularly convenient to work in the ADM formalism |18j since it keeps the invariance under 
spatial diffeomorphisms manifest: only quantities that are covariant under spatial transformations 
appear in the equations. The perturbed four- dimensional metric in ADM variables is 

ds 2 = -N 2 dt 2 + h ij {dx i + N i dt)(dx j + N j dt) , (A.l) 

where N(pc, t) and iVj(x, t) are the 'lapse' and 'shift' functions, respectively, and hy is the induced 

metric on three-dimensional hypersurfaces of constant time t. The geometry of these hypersurface 

(3) 

is characterized by the intrinsic curvature Rh , i.e. the Ricci tensor of the induced 3-metric, and 
the extrinsic curvature Kif 



Eu = NKi 



1 



(A.2) 



where V is the covariant derivative associated with hij. We will use these geometrical quantities 
to define the most general action with unbroken spatial diffeomeophisms. As usual in effective 
field theories, this can be done in a low-energy expansion of the fields and their derivatives: We 

(3) 

notice that N is a scalar with zero derivatives acting on it, E^j is a one-derivative object and R^-' 
contains two derivatives. 

A.l Zero-Derivative Action 

At the zero-derivative level in ADM variables the only quadratic operator is (5N) 2 . Let us 
therefore study the action 



S = I d 4 x^=5 



i? (4) - M 2 i (-^H + 3H 2 +H\+ M\t){5Nf 



(A.3) 



2 

where R^ = R^ + N- 2 (E i ^E ij - Ef). We fix spatial diffeomorph isms by the gauge choice 

hij = a 2 (t)(l + 2()5 ij . (A.4) 

The Einstein constraint equations express and iVj in terms of £. In particular, variation of the 
action with respect to N gives the Hamiltonian constraint 



Ml 



R {3) - A ' 2 (/•."' /-',; - E\ 2 )\ - Mpj (-^H + 3H 2 + 2M 4 5N = , (A. 



2 

while the variation with respect to N l gives the momentum constraint 

1 

N 



5) 



M 2 l V j (-(Ei-5iE k k ))=Q. (A.6) 
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To get the quadratic action we only have to solve these constraints at first order in £. Defining 
N = 1 + JVi and JV* = doji, we find 



H 



2M 4 H \ ■ V 2 C 
M^H 2 ~ IP J C " ' 



(A.7) 
(A.8) 



Substituting this back into the action A.3 we get 



S 2 = Mpj / d 6 xdta\t) 



A(t)( 2 + B(t) 



where 



6 H 
A{t) = -= and B(t) = -e 



H 2 ' 



(A.9) 



(A.10) 



The dynamics of £ is hence determined by two functions of time: the Hubble expansion rate H(t) 
and the 'speed of sound' 

-7T = 1 • (A.ll) 

c 2 M 2 ,H { 1 

Defining a new time variable y via dy = (c s /a)dt, the action becomes 

S 2 = M 2 l fd^xdyq 2 [(C) 2 - (d,C) 2 ] I , where q 2 = — 



(A.12) 



This action is the starting point for much of our discussion in the main text of this paper. 
A. 2 Higher-Derivative Terms 

At next order in the derivative expansion we must include terms that involve the extrinsic cur- 
vature E, L j. We will be interested in fluctuations around the background Em = a 2 H5ij, i.e. 



&E>ij — -^ij ^ij ' 



(A.13) 



There are now three additional quadratic operators: 5N5E\, 5E t: >5Eij and (SE^) 2 . In order to 
simplify the calculations, we will treat the operators one at a time. 



(5E\f Term 



Let us first add the operator (5E\) 2 to the action (A.3), 
S 



I 



d X\f^g 



R {4) - AfJ (J^H + 3H 2 + H^j + M 4 (t){5N) 2 - M 2 (t){5E\f 



(A.14) 



The Hamiltonian constraint is the same as in (|A.5|), while the momentum constraint becomes 

(A.15) 



M 2 x V j ( -{E\ - 5lE k k ) ) - 2M 2 V i E k k = 
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The first order solutions to the constraint equations now are 

M^H ■ ( + M^M 2 ■ V 2 (/a 2 



M 4 if 2 + M 2 M 4 
(-9M2jM 2 i? 2 + MpjM 4 - M 4 ii)C - Af 4 # • V 2 C/a 2 



Substituting this back into the action we get [8] 
S 2 = Mpj / d 3 xdta 3 (t) 



i(t)C 2 + S(t)^ + C(t)^ 



where 



B(t) 
C(t) 



M^M 4 - 9M 2 if 2 - M^H) 
M^H 2 + Af 2 M 4 



(A/ 4 # 2 + M 2 M 4 ) 2 



3M^M 2 i? 4 + AP(M 4 - M^H)(M^M l - Mt x H) 



4ji>2 



H' Z (M 4 M^M 2 + Af^fT) - M^M'^HH 



6 Ji/r2 i 



M^M 2 



M^H 2 + M 2 M 4 ' 



In Appendix [B] we will analyze certain interesting limits of this action. 



(A.16) 
(A.17) 

(A.18) 
(A.19) 

(A.20) 
(A.21) 



5N5E) Term 



Next, we add the operator 5N5E\ to the action (A. 3) 
d 4 Xy/—g 



S 



Af 



1 



- f ^ff + 3fi 2 + H ) + M 4 (i)(<5iV) 2 - M 3 (t)5NSE\ 



Performing a similar analysis as before we find [8] 



S 2 = Mpj / d 3 xd£a 3 (i) 



A(t)C 2 + B(i) 



($0 
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where 



A(t) 



B(t) 



■UM^HM 3 + 4M 4 M 2 1 - 4Af 4 1 ff 



(-M 3 + 2M 2 x Hf 



AT 6 - 2M 2 X HM 3 + 4M 4 if 
(—AT 3 + 2Af 2 if) 2 



(A.22) 



(A.23) 



(A.24) 



(A.25) 



This action is of the same form as (A. 12) but with a different normalization and a sound speed 
given by 



c 2 = 



B Af 6 — 2M 2 1 i?M 3 + 4M 4 ii~ 



(A.26) 



Notice that the action is defined by three functions of time: a(t), M(t) and M(t). 
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B Scale-Invariance from Higher-Derivative Terms? 



Could the actions (A. 18) and (A. 23) allow for scale-invariant fluctuations in non-de Sitter back- 
grounds without leading to a strong coupling problem? The goal of this appendix is to address 
this question. Since a comprehensive analysis of the full range of possibilities is beyond the scope 
of this paper, we will focus on a number of interesting limiting cases. 

5N6E\ Term 

For notational convenience, we define the following dimensionless ratios 

X= " and Y = - \ „ , (B.l) 

M 2 H 2 2M 2 H 2 ' V ' 

pi pi 

where, despite appearance, X and Y can both be negative. In this parameterization we find 

A = (e + X -QY)(1-Yy 2 , (B.2) 
c 2 = [ e + Y -Y 2 )(e + X -6Y)' 1 . (B.3) 

We assume M ~ M and H <C M, so that |X| ^> \ Y\. We will furthermore restrict the discussion 
to the following two limits: 

i) |M 4 | < M^\H\ O \X\ < e , 
it) |M 4 | > M 2 X \H\ ^ \X\ » e . 



In the first case, |Y| < |X| < e, ( |B~2| ) and pO] ) simplify to 



i -> e (1 - F) 2 and -»■ 1 . (B.4) 



For \Y\ <C 1, this reproduces the canonical action of A — > e. In this case, we know that 
only slow-roll inflation leads to a weakly-coupled solution. For |y| > 1, the solution differs from 
slow-roll inflation, but it seems hard to avoid the strong coupling problem from interactions 
proportional toe^>|X|>|Y|>l. Nevertheless, it is possible to find a scale-invariant attractor 
solution by choosing 1 — Y oc r, while keeping e large and time-independent. This would realize a 
contracting ekpyrotic universe. However, in this case metric perturbations 5N oc (1 — Y) _1 oc r^ 1 
become large, and interactions mediated by the mixing with gravity become strongly coupled. 
The more interesting limit therefore is the second case, \X\ 3> e, for which we obtain 

A -> X(l-Y)~ 2 and c 2 s -> ^(e + Y -Y 2 ) . (B.5) 

X 

To avoid a ghost instability we require A > and hence X > 0. For | Y| < min{l, e}, the extrinsic 
curvature corrections proportional to Y are subdominant and the theory reduces to 

c 2 -> 4 < 1 and i -»■ X = 4 • ( B -6) 
X c 2 

This is the case that we studied in the main text and which produced inflation as the only 
weakly-coupled solution. 
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The next parameter region to consider isl<|y|<^e<^|X|. The coefficients of the quadratic 
action become 



~ J«l and A -* J^-^-Y)- 2 - (B.7) 



In this limit, interactions proportional to X have to be controlled to avoid a strong coupling 
problem. For instance, a particularly dangerous interaction is 



C 



«^) 2 !L(i_y)2 C . (Bi 



c s 

We require Y to be sufficiently different from 1, in order to avoid the strong coupling problem 
discussed below Eqn. (B.4). Since e is bounded from below, e>l, and c s is bounded from above, 
c s < 1, we cannot afford much time- variation of both e and c s while keeping the fluctuations 
weakly coupled. We therefore consider e ~ const and c s ~ const. The condition for scale 
invariance then implies 

a 2 A c s = ^ (1 - Y) ~ 2 oc y- 2 oc r~ 2 . (B.9) 

C s 

For constant e, we get a power law solution for a(r), 



2 - => aoc r 1/(e - 1] , (B.10) 



and hence 

1-7oct h A~ t . (B.ll) 



This quickly leads to a strong coupling problem as Y approaches 1 — see the discussion below (B.4). 



The last case we need to discuss is: 1 > F > e ( 



1.2 



Y * „ Y 

X 



c 2 . -»• ^7 < 1 and I -> X = ^ . (B.12) 



As before, interactions proportional to c s 2 restrict us to solutions with c s ~ const. Scale invari- 
ance then requires 

a 2 Y M 4 1 

a 2 Ac s ~ — ocr- =► Y oc X = ^ oc ^ , (B.13) 

For ar = const, this reproduces the near-de Sitter solutions of II 2| . For a ~ const., Y 
grows rapidly and we again have to worry about a strong coupling problem as Y approaches 1. 
Moreover, the time-derivative of (B.13) implies 

/dlnax- 1 

2 « S3 - 1 -(^) ■ (R14 > 



12 We disregard the case \Y\ > 1 because it leads to c 2 . < 0, signaling that this limit is generically plagued by 
gradient instabilities. These instabilities may be cured by adding the (8E\) 2 operators, but only at the expense of 
having the two higher-derivative terms be equally important near the horizon scale, as described in [12! . While this 
is quite a reasonable and technically justified assumption in the case where all the parameters are approximately 
time-independent as in inflation [12] . we do not consider this an interesting scenario when all of the parameters 
are strongly time-dependent. 
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where we used that e <C 1. For |ejy| -C 1, (B.14) describes perturbative deviations from the 



near-de Sitter solutions of 112] . Finally, we consider the limit |ej\/| 3> 1, for which (B.14) 
becomes 

/dlna\ _1 

2eM~-(-77-T) • (B-15) 



Vcilnr 

We have to worry that the time-dependence of M(t) induces a strong coupling problem via 
interactions of the forrrF^I 



CC 2 



(B.16) 



We see that cm cannot be too large. If it starts large and increases with time, we generate large 
non-Gaussianities. If it decreases with time, it quickly enters the regime \cm\ ^ 1- Hence, just 
as in the case of large e, we impose €m — const, in order to keep (B.16) under control, while 



remaining in the non- inflationary regime \cm\ 3> 1- In this case we can integrate (B.15) to get 

-l/2e M 



However, this power law solution for the scale factor implies e 
inconsistent with our original assumption e <C 1. 



(B.17) 

2 1 cm | 3> 1) which is of course 



Summary. Although we haven't given a completely comprehensive analysis of the action ( A. 23 ), 



our analysis of the most promising limiting cases has only revealed the near-de Sitter models 
of [3 H2] as backgrounds with weakly-coupled scale-invariant fluctuations. 



(6E 



i\2 



Term 



To analyze (A. 18), it is convenient to define 

x M ' 



and 



Y 



M 2 



(B.18) 



B 

C 



(l + XY)- 2 
-(1 + XY)- 1 



+ Yrje 



(B.19) 
(B.20) 

(B.21) 



so that the coefficients of the action can be written as 
A = {1 + XY)- 1 (X - 9Y + e) , 

-3Y + XY[XY + e(l + Y)] + [XY 
Y 

H 2 ' 

As before, we assume M ~ M and H/M <C 1, so that \X\ 3> Moreover, we impose \Y\ < 1 
because the cutoff of the effective theory shouldn't exceed the Planck scale. We notice that the 
action reduces to the correct limit of slow-roll inflation — A — > e, = —B A -1 — > 1 and C — > — 
for l3>e^>|X|>|y|. In order to be maximally different from slow-roll inflation we will restrict 
our discussion to the limit \X\ 3> e: 



A 
B 

C 



[l + XY)' 1 X 
'l + XY)- 2 

-(1 + XY)- 1 



-3Y + XY[XY + e] + [XY 
Y 

H 2 ' 



+ Yrje 



(B.22) 
(B.23) 

(B.24) 



13 Here we derive the cubic interaction from the decoupled 7r-lagrangian [9] ignoring the mixing with gravity. We 
expect this to give a lower bound on the level of strong coupling. 



23 



For Y = this reproduces models with small speed of sound — A — > X, c 2 



x 



< 1 and C 0. 



For <C 1, and assuming the higher-derivative term to dominate at horizon crossing (this 

has to be checked a posteriori for each solution), the action becomes 



where 



S 2 = M^J d^xdrq 
0/ . a 2 M 4 



(C 



l\2 



M 2 tf 2 



and /u (r) = 



a 2 M 4 
M 2 ' 



(B.25) 



(B.26) 



We require a slightly new treatment to determine the conditions under which (B.25) leads to 
scale-invariant two-point correlations: We first consider the equation of motion of the canonically- 
normalized field v = v2 q( 



v'l + W 



.2 Q" 



Vk = , 



Modes freeze when 



or 



where 



k 4 

oj 2 (k,r) = . 

I 1 



(B.27) 



(B.28) 



At early times, i.e. when all modes are deep inside the horizon, the equation of motion reduces 
to 



v'l + uj 2 v k = 



or 



d 

7h 



1 



UJ 



UJ 



UJ \v k \ 



Averaging over many oscillations, (|f^.| 2 ) 



(u 2 \v k \ 



1 dE _ 1 doj 
E dr uj dr 



2E, we find 
E 



or 



const. 



UJ 



This implies 



ICI 2 oc 



*L0 



l {q"/q) 1 / 2 



Assuming q to be a power law in conformal time, 



q oc r 



this becomes 



ICI 2 oc 



1 



T 



2n-l 



where 



k. oc 



If we parameterize the time-dependence of u as a power law, 

rP 



we find k 4 oc r 2<yP ^ and hence 



1 — 4 "- 2 

|C| oc k p- 1 



(B.29) 

(B.30) 
(B.31) 

(B.32) 
(B.33) 

(B.34) 
(B.35) 



We see that scale-invariance, |£| 2 oc k 3 , requires the following algebraic relation between the 
power law indices 

An = 3p - 1 . (B.36) 
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The time-dependence of the canonically-normalized field on superhorizon scales is v oc r ra and 
r~ n+1 . Hence, the two solutions for £ ~ vq _1 scale as r° and r _2n+1 . For n < g the growing 
mode of £ will be a constant and the background is an attractor. We will restrict to that case. 
Let us further simplify the treatment by considering solutions with constant e and hence 

4n + l 



(x r 1 /^ x ). The conditions for scale-invariance, q oc r n and /u oc t 3 , then imply 



M 2 M 3 1 
m oc — oc a 



2e M + 3e 



M 



-5e . 



(B.37) 



We have therefore found a whole family of solutions with scale-invariant two-point functions. This 
includes ghost inflation [TU]: e = and em = = 0, i.e. de Sitter backgrounds, a = —1/(Ht), 
with M = M = const. We also find de Sitter backgrounds with possibly large time-variations 
in the coefficients of the effective action: e = and cm = —\^m- Finally, the solutions (B.37) 
contain non-inflationary cases with e = const. ^> 1 and {|cm|) \ € m\} = const. 3> 1. 

Let us look at the last solution in a bit more detail. For \X\ 3> e = const. > 1 > \Y\ and 
I -XT' I <C 1, the action (A. 23) simplifies considerably 

[(O a -3(*C) 2 

where 



S 2 = Mpj / d 3 xdrg 



2 (d 2 C) 2 



(B.38) 



e 

X 



a 2 X 



a 2 e 



(aH)' 



X 
Y 



The /c 4 -term dominates at horizon crossing if 



22 k t 



where 



n[n 



V 2 



ran 



1) e 2 (a#) z 



|X| 

IFF 



Substituting ( |B.39[ ) and ( |B.4ip into QB.40[ ) we find 

1 > ixyi > 



Since we restricted to n < ^ for attractor solutions, this implies 

n< -(1 -V5) ~ -0.6 . 



(B.39) 



(B.40) 



(B.41) 



(B.42) 



(B.43) 



This forces the coefficients of the action to be strongly time-dependent. We should therefore 
check the time-dependence of interactions. Consider, for example, 



and compare it to 



£ ('(a,() 2 



-(C) 2 



2 2 r/-i\2 



MW (C 



(B.44) 
(B.45) 
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Using oj = k 2 1 u, we find 

W_/. :( = !^! (j (B . 46) 



£ (c , )2 aH^ |y|V2 



where 

v v , _ . _ 



XI IXL; /r\f(n+l) 



(B.47) 



Here, we have used that u oc r 3 + and a-ff oc r 1 . Since IXL 3> \Yh, this is a large interac- 



tion at some fiducial time r%. Moreover, it grows, unless n > — 1. Together with (B.43), this 
essentially rules out the scenario. This case is therefore another example where two-point cor- 
relations are perfectly scale-invariant, but higher-order correlations are forced to have a strong 
time-dependence. The fluctuations therefore again can't be weakly coupled over a sufficiently 
large range of scales. 

Finally, we consider the opposite limit \XY\ S> e ~ const. > 1. In this case the action is of 



the same form as (B.38), but the coefficients now being 

c^-Y, q 2 ^y = -^, f^(aHf*. (B.48) 
Y Cg Y 

The /c 4 -term dominates at horizon crossing if 

e < \XY\ < e 2 \n(n- 1)| . (B.49) 

For e > 1 this is easily satisfied. Note that our previous conditions for scale-invariance, q oc r n 

and n oc r 3 , still apply. This time an interaction proportional to M creates the most obvious 
problems: 



c 



(CO 



A2 



-M 



XY(, (B.50) 



T\f(2-n) 



\XY\ = \XY\i ( — 1 3 ' . (B.51) 



where 

\XY\ = \XY\t 

In order for this large interaction to be time-independent we require n ~ 2 (for n > 2 the 
theory quickly becomes strongly coupled, whereas for n < 2 we quickly exit the regime of interest 
\XY\ 3> e.). However, this is larger than our upper bound on attractor solutions n < \. We 
therefore don't find a successful attractor solution in the regime \XY\ 3> e. 

Summary. We have reproduced ghost inflation, but found no successful non-de Sitter back- 
grounds with scale-invariant weakly-coupled fluctuations. 
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C Comments on Non-Attractor Solutions 



In this appendix we consider the second type of backgrounds that lead to scale-invariant fluctu- 
ations: 



Since the curvature perturbations £ in this case evolve outside of the horizon as £ oc y , we can't 
make definitive predictions for observations without specifying the evolution of £ after the scaling 
phase. We will work throughout with the simplifying assumption that £ freezes immediately after 
the scale-invariant phase ends. The history of the FRW background during the scaling regime 
is then sufficient to predict the primordial power spectrum relevant for observations. Many of 
these results have previously been derived by Khoury and Piazza [10] . This section is a concise 
summary of their analysis with a few corrections and clarifications. 

C.l Basic Challenges for a Predictive Theory 

The non-attractor nature of the background implies that additional physics has to be added 
before and after the scaling phase. This significantly complicates even the minimal scenario. We 
will first indicate the basic model-building challenges, before adding more technical details and 
an analysis of the strong coupling problem in the following sections. 

The basic timeline of the universe is: 

• Pre-scaling phase 

Since inhomogeneities grow exponentially during the scaling regime they need to be expo- 
nentially small at the beginning of the scaling regime. To avoid an enormous fine-tuning 
of the initial conditions requires a pre-scaling phase with the right physical characteris- 
tics to set the initial conditions dynamically, i.e. the universe needs to be cleaned from 
pre-existing inhomogeneities to extremely high precision. This is similar to the cleaning 
mechanism proposed for New Ekpyrotic Cosmology [17] . 

• Scaling phase 

The background then develops a specific time-evolution of the speed of sound that allows 
for the generation of scale-invariant curvature perturbations. If a(y) and c s (y) are power 
law solutions, scale-invariance can only be achieved if c s changes rapidly with y. This is 
consistent with the result that curvature perturbations in a contracting spacetime have 
a very blue spectrum if c s = 1 [T3j. Just as in the case of the exponentially changing e 
parameter above, the exponential change of c s is reason to worry about the strong coupling 
problem. However, now there are two independent functions of time - a(y) and c s (y) - 
as well as the time-dependent curvature perturbation so we expect more freedom to 
evade constraints. 

• Post-scaling attractor 

A crucial challenge for the non-attractor cases is to make contact with late-time cosmo- 
logical observables. The time-evolution of curvature fluctuations needs to be followed even 



2 
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after the modes exit the horizon. During the scaling phase this time-evolution is fixed, 
however, after the scaling phase there is some freedom in matching to an attractor solution 
that connects continuously to the FRW expansion after the bounce. To avoid introducing 
an element of unpredictivity, we follow Khoury and Piazza [10] in assuming that the back- 
ground transitions instantaneously from the scaling phase to a post-scaling attractor, i.e. £ 
(and c s ) freeze immediately after the scaling phase. Observables then only depend on the 
physics during the scaling regime. 

• Bounce 

As usual in contracting models, the background has to bounce to connect to the expanding 
FRW phase. If the background has approached the post-scaling attractor before the bounce 
then the scale-invariance in the pre-bounce fluctuations will be preserved in the post-bounce 
fluctuations [14J (unless the unknown physics of the bounce is sensitive to exponentially 
small decaying modes). 

Given these minimal model-building requirements, we next discuss the basic constraints aris- 
ing from scale-invariance and weak coupling. 



C.2 Background Solution 

To characterize the time-dependence of the background during the scaling regime, we define the 
following two parameters 

< — - and e s 



He, 



In terms of derivatives with respect to y (d/dy =') these can be expressed as 



(lnc s )' 



and 



2 - e, - e = 



a" la 



(In a)' " ^ s " (a'/a) 2 ' 

Assuming that e and e s are constant gives power law solutions [10] 

a(y) ~ {-yf and c s (y) ~ {-y) pts , with p 



e s + e-l 



and 



~ (-y) 6 , with 5 = 



2-e s 
e s + e-l 



Comparing (C.5) to (C.l) we find 

Case II 

and 



5 = 4 



e s = -(3-2e) 
5 



5 S — 4. 

i(y) ~ {~y) 1+e and c s (y) ~ (-y) 1+* 



(C.2) 
(C.3) 

(C.4) 
(C.5) 

(C.6) 
(C.7) 



Via the algebraic relation between e and e s in (C.6) the combined time-dependence of a and c s 
ensures scale-invariance of the two-point function. In the analysis of the strong coupling problem 
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the time-dependence of the term C/ C s between horizon exit and the end of the scaling phase will 
be essential, 



c 


Cend 


( y y 




~~ C 2 

s,end 





with 



P = 



5(6-3) 
1 + e 



(C.8) 



where Cend and c S)en d are constrained by observations: Cend ~ 10 5 and c Sjen d > 0.01. 



Digression: spacetime during the scaling phase 

Since inhomogeneities grow during the non-attractor phase it is interesting to determine what the 
spacetime would look like. To first order in £, we find the following metric |14j 



ds 2 -- 

The intrinsic curvature is 



1 + 2 — 


dt 2 + di 


. H. 





c? V 2 ^ H 



dtdx 1 + a 2 (l + 2C)% dx l dx 3 . 



H 2 ~ (aH) 2 ' 



In the background (C.7) this has the following time-dependence 



H 2 



(kyy\v 



(C.9) 



(CIO) 



(C.ll) 



The intrinsic curvature perturbation is hence suppressed on superhorizon scales \ky\ <C 1. Furthermore, 
unless ft is very negative it decreases with time. From the metric we can also compute the extrinsic 
curvature 



1 



K\j = -g lk {g kl - d k g 0j - d j9m ) 



To first order in Q, we find 



5K- 



3 _ C 5 i djdj 



H H 



lL 1 (3)r 



(C.12) 
(C.13) 



where K % ^ — HS* is the extrinsic curvature of the background. Since 1 has the same time-dependence 

as C, 

H (may 

the isotropic component of the extrinsic curvature tensor grows exponentially relative to the background. 
The time-dependence of the anisotropic component of the extrinsic curvature follows from 



IC_ 

c 2 H 



\y 







(C.15) 



Only for /3 not too negative will the anisotropic extrinsic curvature be under control. 
C.3 Strong Coupling Constraints 

When is the time-evolution of a and c s consistent with weak coupling over a large range of scales? 
As before, we require that the theory is weakly coupled at horizon crossing, i.e. 



X 



0({l,e,W) 



C 



< 1 • 



(C.16) 
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This guarantees that quantum mechanical loop effects are small and the perturbative computation 
of the primordial fluctuations is under control. Furthermore, the growth of £ on superhorizon 
scales can lead to the classical growth of non-linearities. Keeping these non-Gaussianities below 
current observational constraints imposes further constraints on the parameter space. 



Quantum mechanical constraints. Eqn. ( C.8 ) describes the size of C%/ C2 at horizon crossing. 



We distinguish the three cases (3 = 0, (3 < and (3 > 0: 
P = 

For e = 3 (i.e. pure kinetic-dominated contraction) the time-dependence of the speed of sound, 
c 2 s oc y~ 3 , exactly compensates for the time-dependence of £ - i.e. both the speed of sound and 
the curvature perturbation are exponentially small at the time when the first modes exit the 
horizon 14 c 2 (yi) < (e AAr )~ 3 ~ e~ 3AJV and £ ~ Ccnde _3AAr , where AiV denotes the number of 
e-folds of contraction between %n and y e nd = e~ AN Vi- The time-independence of Q/c 2 s ~ 10~ 5 for 
f3 = implies the absence of the quantum mechanical strong coupling problem. However, below 
we show that the superhorizon generation of non-Gaussianities in this case leads to a classical 
strong coupling problem that rules out the j3 ~ limit of the parameter space. 

P < 

For e < 3 the fluctuations are more weakly coupled at the initial time yi than at the final time 
y en d- Hence, if the weak coupling condition is satisfied for the mode that exits the horizon at 
Vend ( as required by observations), then it will automatically be satisfied by all modes that exit 
the horizon at earlier times y < y en d- 

13 > 

For e > 3 the fluctuations are more strongly coupled at the initial time yi than at the final time 
Vend- This fact implies an upper limit on e arising from the fact that we want to avoid a strong 
coupling problem for the largest CMB scales, which cross the horizon at yi. Quantitatively, the 
constraint 

X( Vi ) = 0(1) • 4^ e^ AN « 1 , (C.17) 
c 

s,end 

implies 

If) 8 

Observations imply at least AiV ~ 10 e-folds of scale-invariant fluctuations. The constraint in 



(C.18) then becomes e < 3.8. However, below we will see that for any value of (3 > (e > 3) the 



non-Gaussianities developed by superhorizon evolution quickly lead to a strong coupling problem. 

Classical constraints. The non-linear growth of £ on superhorizon scales leads to an additional 
accumulation of non-Gaussianities. In this subsection we identify the regime of parameter space 



14 The model therefore features a couple of very small (but technically natural) numbers: e.g. the sound speed at 
the beginning of the scaling phase is (? s (yi) = 10~ 13 for AiV = 10 and <? s (yi) = 10 -65 for AiV = 50. Furthermore, 
using H^cJ 1 ~ y~ 6 we find H(jji) ~ e~ 10 ~ ~r AN M p i, i.e. the Hubble parameter at the beginning of the scaling 
regime is H( yi ) = 10~ 21 Af p i for AiV = 10 and H{ Vl ) = 10 _86 M p i for AiV = 50. 



30 



for which those non-Gaussianities remain small enough to be consistent with observations. In 



Appendix |D| we compute the exact bispectra for the background of {C.2 following Khoury and 
Piazza ITUl: 



(Ck 1 Ck 2 Ck 3 ) = (2vr) 7 5(k 1 +k 2 + k 3 



EI A 3 



A . 



We will use two different measures for the amount of non-Gaussianity: 

f - 40 1 V A 

6 K a fci— »0 

for non-Gaussianity in the squeezed limit, and 

/nl = 30 — g Ak 1= k 2 =k 3 , 



(C.19) 



(C.20) 



(C.21) 



for non-Gaussianity of general shape. Here, we have defined K = k\ + k<i + . As before, we will 
describe our results in term of deviations from the sweet spot e = 3, or (3 = 0. The details are 
presented in Appendix [Dl Here we summarize our findings: 



(3 = 

In the limit (3 = simple analytical expressions for the non-Gaussianities can be obtained. For 
instance, the interaction Q J produces the following non-Gaussianity in the squeezed limit 



c 3 

lim f^, 

/3^o NL 
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(fP - 1) • log(k/k end ) + 



"s end 



(C.22) 



where is an order unity coefficient related to coefficient of the (<5g 00 ) 3 term in the effective 



D). Unless one "tunes' 



Lagrangian of [9] (see Appendix 

non-Gaussianity. This limits the scaling regime to a 

" A7V , Eqn. \C.22\ becomes 



15 



fx^ this gives a large scale-dependent 
::ew e-folds. Specifically, for c sen d ~ 1, 



1| ~ O(l), and fe/feend = e 



The observational limit |/nl| i$ 100 then implies 



24 AN . 



(C.23) 



AN < 5 . 



(C.24) 



This is of course inconsistent with the observationally required range of scale-invariant modes, 
AA^ > 10. This result is not a quantum disaster as before (in the sense of > 1 at horizon 
crossing), but only an observational fact about our universe. Moreover, we should remark that this 
constraint can be avoided if \ — 1| <^ 1. In fact, if the background is sourced by a Dirac-Born- 



Infeld (DBI) action then fjp = 1 and the dangerous contribution in (C.22) vanishes. In that case, 



the leading non-Gaussianity comes from two non-local interactions (see Appendix|D[). Individually 



'This "tuning" may have a microphysical justification. In fact, curiou sly, if the background arises from a 



Dirac-Born-Infeld (DBI) action then /j? = 1 and the contribution in (C.22 1 vanishes 
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each of these interactions produces a large scale-dependent non-Gaussianity. However, in the 
squeezed limit the two terms cancel precisely up to terms proportional to f3, 

lim + /£) = 0((3) . (C.25) 

Remarkably, for (3 = the combined non-Gaussianity produced by these interactions is therefore 
identically zero in the squeezed limit. However, this cancelation doesn't hold in the equilateral 
configuration, as we see from the following result: 

lim + /£) = -J-^log(fc/fc end ) + 4^- « 10AJV. (C.26) 

''s end *\s end 

At best this is marginally consistent with the observational limit /nl ^ 300 |12j . 

The above proves that the special case (3 = (corresponding to pure kinetic energy domi- 
nation), although consistent with £3 <C £2 at horizon crossing, is ruled out by the size of the 
non-Gaussianity arising from the superhorizon evolution of £ (unless fy' is close to 1.2 — see 
below) . 



/nl 








Figure 4: Plot of /nl vs. /? for j/j = e 10 y C nd (blue) and yi = e 5t) y cn d (red). We show two cases: f x = 1 
(DBI) (=*> /nl > 0) and (/^ — 1) = 1 (=>■ /nl < 0/ The piot implies an upper limit on (3 but no lower 
limit. This difference from the result of f J Of arises from corrections to their result for the bispectrum (see 
Appendix ITj|). 



(3 < 

For general (3 we present the result graphically in Figure [4] (for specific values of fjP)- The 
plot confirms that the magnitude of /nl is inconsistent with observational limits for (3 > 0. 
However, it also shows that /nl(/3) is strongly suppressed for negative (3 (or e < 3). For generic 



32 



f^' — i.e. \ f^ — 1| = 1 — we find a new upper limit for the equation of state during the scaling 
regime: 

< -0.3 =>- e < 2.77 . (C.27) 

For e < 2.77 the non-Gaussianity is still consistent with present observational constraints. More- 
over, if fx is tuned to lie between 1.1 and 1.3, the non-Gaussianity is sufficiently small for all 
negative j3. 

It is interesting to compare this upper bound on e with a lower bound obtained by requiring 
that the background during the scaling phase solves the flatness, horizon and isotropy problems 
of standard big bang cosmology. Consider the Friedmann equation 

H --^ + 3Mjv^ + ^ + ^J + 3A5^' (C - 28) 

where we have allowed for contributions from curvature, matter, radiation and anisotropics, in 
addition to the fluid X which sources the background during the scaling regime. We see that e > 3 
is required such that pre-existing anisotropics don't come to dominate the universe. Interestingly, 



this is inconsistent with the bound in (C.27). This implies that in order to make the scenario 
viable, a phase has to be added before the scaling regime which reduces anisotropics to such a 
low level that their subsequent growth remains harmless. This is to be contrasted with Case I 
(inflation) where backgrounds with e < 1 automatically solve the flatness problem without added 
features. 



D Bispectra for Non-Attractor Solutions 



In this appendix we compute the complete bispectra of curvature fluctuations in the non-attractor 
solutions of Appendix [C| following the methodology developed in [TOj . Our starting point is the 
cubic action for £: 
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+ 



+ 



dtd 3 x 



2H 2 e 



3 ci 



fx -(I- c 



a 3 e 



ac 



(e - 2e s + 1 



a 3 e d 



C(ftC) 2 -2a^c(ftC)(5a) 



2c? dt 



^ ( 4 ) C 2 C + ~{diO{d iX )d 2 x + f (d 2 C)(d tX ) 2 + 2/(0 6Cl 



2a 



4a 



(D.l) 



where dots denote derivatives with respect to proper time t, di is a spatial derivative, and x 1S 
defined as 



a 2 e 



C- 



(D.2) 



16 This action follows from the effective Lagrangian of Ref. (9 by including the effects of mixing with metric 
perturbations [15]. The order-one parameter fx is model-dependent — in the effective theory of inflation it arises 
from the coupling of the (5g 00 ) 3 term of the action — (see [51 [TO) for the precise definition). If the action (D.l l is 
derived from the DB1 action then /° BI = (1 - c 2 s ) [TDj or C3 131 = |(1 - c 2 3 ) 12 ]. Here and in the following we have 
set M p i = 1. 
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£ f 

The term denotes the variation of the quadratic action with respect to the perturbation £: 



6C 2 



(D.3) 



Finally, we have defined the function 

/(C) 

1 



1 



4cT ^ c 2 H^ ^ Aa 2 H 2[ 



+ d- 2 (d i d j (d i (d j ())] 



2a 2 H 



[($C)($x)-3-Wi($C0,-x))] 



(D.4) 



where 5 2 is the inverse Laplacian. 

In an impressive tour the force Khoury and Piazza |10j computed the bispectra for the action 



(D.l ). In this appendix we reproduce their results (with small corrections) for the case of constant 
e, e s and 63. We will use the standard expression for the three-point function in the in-in 
formalism: 

(D.5) 
(D.6) 



(C kl Ck 2 Ck 3 ) = -i / dy / ([C(0,k 1 )C(0,k 2 )C(0,k 3 ),^ in t(y / )]) • 
J —00 

The operators £ are expanded in creation and annihilation operators in the usual way: 



k ' 



where [ak, S_ k/ ] = (2w) 6(k + k') and the mode functions are [10] 
Finally, we recall that the background satisfies 



H 

'Jl 2 



H, 



end / V 



s end 



and 



H 

~Yf2 ~ 1/2 



H en d ( y 



(D.7) 



CD.J 



where /3 = 5(e - 3)(e + l) -1 . 



D.l Shapes of Non-Gaussianity 

With these ingredients we compute the bispectra associated with the individual interactions in 



(D.l). Defining K = k\ + hi + k% and 



we find: 



(C kl Ck 2 Ck 3 > = (2^) 7 <5(k 1 + k 2 + k 3 



_3_ 



A, 



(D.9) 



cos 



TTi3 



.(II) _ + /.(o) " 
^C 3 lOc 2 __A Jx J)(K\y ead \)P 



r(/3 + 3)M|M + 3r(/3 + 2) ^3 ^ ^. 



9K 3 



(/3-2)(/3-3) 



cos 



(K\y end \r 



/3 2 - 2/3 + 3 3 



p-1 



+ 



/3 
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where we defined fx = + /x^ c s> 



A 



(II) 

CC 2 



+ 



(e-3) 
4c Lnd 



cos 



7T/3 



{K\ VwA \)P 



r(/3 + i) 



2 + /3 



e ^1 



+ 3r(/3) [ 2 fciA^ + 9fcifc 2 fc 3 + j^z 



9K 3 



(/3-2)(/3-3) 



+ 3T((3-1)K ^2K 2 + 5^2 hk^j 



3 l 2K^k t k J+ ^k lk ]-K^^- 3 

i<j v£j 



(/3-2)(/3-3) 



and 



.4 



(II) 



13 e - 7 cos 



7T/3 



F r(/3 + i) K>f 



The following two terms are from the non-local interactions in Si n t: 



Kj 



A (II) 

c,dcd x 



4 c 2 

^ c s end 



cos^r(/3 + l) 

(^|y c „d|)^ 



2Kk 1 k 2 k 3 



^kfk^ kfk 3 



+ 3 



'cos^r(/3) i N 
C^dF"^ 



9hk 2 k 3 -2^^+ 



i¥=3 



/ cosfr(/3-l) 1 



" 2 E fc I - 9 E fci/c J + 



1 2 3 \ 



E fc2fc i E fc * 4 
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cosf T(/3) + 0-1 



(^bendD^-2 /3(/3-2) 3 



if 3 

0-3 
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A™ = 



32^ end 



cos^r(/3 + l) 



{(3 + 4) \J2 kj ~ E kik ) + 2felfc2fe3 



i2 kl k 2 k 3 - A e ^ + ^^(E % 5 - 2 E % 4 + E k > k l 



+ 



i+3 



-54hk 2 k 3 - 15 E k i k j - 6 E fc i + 



k\k 2 k 3 



3E^ 



3 e *w + 2 E - 4 E k n + p^p E + E W - 2 E fc4fc 

i^j i^j i<j / 1 2 3 V i^j i^j i^j 



( cos ^r(/3 - 1) 1 
3K [ ,Z „ 7 + 



-6 e k ) - 18 E fci/c 



' hk 2 k 3 



(K\y end \)? (3j 
E k j ) + k\k\k\ I 2 ^ ~ 4 ^ + 3 fc ^J ~ 3 E 



E + E ^ 



cos 



7T/3 



r(/3) + /? - i i 



{K\y enA \fP-2 0(13-2) 3 



if 3 
0-3 



h 2 h 2 h 2 



E fc i E fc2 ^i 



Finally, there is a term associated with the field redefinition of £: 



(„) _ 3(1 + 6) 



redef 



(D.10) 



We will use two different measures for the amount of non-Gaussianity: 

/NL = "TT "TT-jT , lim n A 

3 iv d fc^o 

for non-Gaussianity in the squeezed limit, and 



/nl = 30 Afc 1= fc 2= fc 3 , 



(D.ll) 



(D.12) 



for non-Gaussianity of general shapes. 



D.2 Squeezed Limit 

We consider the squeezed limit (k\ — > 0, k 2 — > k 3 = y) of the above bispectra. For the first 
interaction we find: 



/nl 



cos ^ T(/3 + 1) 



"s end 
+ 



(K\y cnd \r 4 



3 9 
+ 



4 /3-1 



G8-2)(/3-3) 



cos 



7T/3 



(^lyendl)' 3 



r(/3) 



/3 2 - 2/3 + 3 3 



/3-1 



+ 



P 



36 



In the limit ft — V this gives 



^3 



lim /nl 



f (o) 



r 2 

Si end 



24 • log(ify e nd) - 19 + 24 7i; + 0(/3) . 



(D.13) 



For generic values of fP ~ 0(1) this gives a relatively large scale-dependent non-Gaussianity. 
However, for the DBI action = 1 and the leading term doesn't contribute. The next interac- 
tion has the following squeezed limit 



f C(?) 2 

J NL 



2 



5(e — 3) r cos 



r(/? + i) + i2r(/3) 



ft 2 - 5ft + 12 
(/3-2)(/3-3) 



+ 78r(/3-l) 



6 2/3 2 - 10/3- 15^ 5 



ft (ft-2)(ft-3) 1^2 



which for small ft becomes 



So /nl 



23 



2 c 2 



Next, we have 



f C(VC) 2 

•/NL 



13e - 7 cos 



s end 



TTii 



+ 0(/3). 



24c^ 



end 



T(ft + 1 



ft + 3 
'/3-1 



which for small ft is 



(D.14) 



(D.15) 



(D.16) 



lim f C( ^ )2 

For the two non-local interactions we find 



c 2 

s end 



+ <W 



J Nh 



30e 



end 



i / cos "frjft) i 

8 \ (K\y end \)P ~ft 
3 



+ 



+ 



cos 



-6 



'cos^r(/3-l) i v 

(^lyendl)' 3 + ft / 

T(ft) + ft-l 1 



/3-3 I (A"|y end |)/ 3 /3-2 /?(/?- 2) 



Jnl 



6 #CdC£>x 
3 



JNL 



It is interesting to consider the limit ft — > (e — > 3). In this case, we get, for example, 



So /nl 



4c 



15 135 

_ (_ 4 + 9 7E ) + _ log(Ky cnd ) + O(ft) . 



(D.17) 



(D.18) 
(D.19) 

(D.20) 



',s end s end 

This looks like a dangerously large amount of non-Gaussianity arising from the superhorizon 



evolution of £. However, from (D.19) we get an equal contribution with opposite sign for ft = 0, 
so that the sum is proportional to ft, 



ft 



c 2 

s end 



(-4 + 97£/) + 91og(ify, 



end; 



0(ft 2 



(D.21) 
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For (3 -C (log(ETy en d)) _1 and c scn d ~ 1 the non-local interactions are therefore less dangerous 
than they appear at first sightj^] However, it is worth pointing out that this cancellation only 
holds in the squeezed limit. For the equilateral configuration, for example, we find 



+ JnlJ 



& . 32 

s end 



log(i^ end ) + -2— + 0{fi) 



end 



(D.22) 



Finally, from the field redefinition we get 



f redef 

Jnl 



1 + e 

' r 2 
u s end 



s end 



1 + 







(D.23) 



In { C.3 we present results for /nl and /nl as a function of /3. 



This cancellation didn't happen in [10] due to algebraic errors in their bispectrum computation. 
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